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Abstract 

We consider pair production in a space-time-dependent background field and 
derive a source term, i.e., production rate in the one-particle phase space. 

Such a source term is required in Boltzmann-equation-based models of quark- 
gluon plasma formation and evolution. We compare the source term derived 
here with the one that has been used in the literature so far. Significant 
differences are observed. 
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I. INTRODUCTION 


In some models of ultrarelativistic nucleus-nucleus collisions, the formation of quark- 
gluon plasma (QGP) is assumed to proceed via creation of a strong colour-electric field 
in the region between the two nuclei receding from each other after the collision, and the 
subsequent decay of the held by parton-pair production according to the well known nonper- 
turbative formula given by Schwinger [^. We shall call these models |^|-0 background-held 
models. They make use of the framework of the relativistic transport theory to study the 
formation and preequilibrium evolution of QGP, with Schwinger’s formula providing a basis 
for derivation of the source term (i.e., production rate in the one-particle phase space) in 
the Boltzmann transport equation. 

In analogy with an ordinary parallel-plates capacitor, the above models have also been 
called colour-plates models. The colour held, however, need not owe its existence entirely 
to the colour plates. A high-density many-body system is often conveniently described in 
terms of a self-consistent mean held and residual short-range interactions. For example, 
Blaizot and lancu P| have considered gauge theories at high temperatures and shown that 
long-wavelength excitations of QGP can be described as collective oscillations of gauge and 
fermionic mean helds. They have presented non-abelian solutions to the equations of motion, 
which correspond to spatially uniform colour oscillations. Importance of the mean held has 
also been stressed in the context of hot and dense hadron gas models of relativistic heavy 
ion collisions; see for example |P]. Whatever may be the origin of the colour held, once it 
is there, it can produce pairs of partons, add to the entropy of the system and ahect its 
subsequent evolution. 

A few remarks on the source term mentioned above are in order. The exact one-loop 
nonperturbative result (Fig. la) that Schwinger derived for pair production in a constant, 
uniform background electric held {E) is given by 


dW 
d^x dt 


aE^ ~ 1 

— exp 


7r2 


mrm 


\eE\ 


( 1 . 1 ) 
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where a = is the hne-structure constant and m is the particle mass. It is instrnctive 


to nnderstand this resnlt in terms of a semiclassical tnnneling across the mass gap This 
resnlt, strictly speaking, is applicable to the parton pair prodnction only in the so-called 
abelian approximation. Secondly, Eq. (1.1) gives the probability per nnit volnme and per 
nnit time, that a pair is prodnced, while the transport eqnation reqnires the sonrce term 
in one-particle phase space. Finally, and more importantly, nnlike the case considered in 
IK] where the constant electric held is maintained by external sonrces, the pairs prodnced 
in the backgronnd-held models derive their energy and momentnm from the held, the held 
becomes time dependent and Eq. (1.1) stops being applicable. 

The resolntion of these problems has so far involved rewriting the right-hand side of Eq. 
(1.1) as an integral over p^ or rri^ = p^ + px being the parton momentnm transverse 
to the direction of the collision axis. The expression so obtained has often been nsed in the 


QGP transport eqnation and it reads 

dW 


d^x dt d^pT 


eE In 


1 — exp 


vrmh 

'W\ 


( 1 , 2 ) 


While the above extension may be justihed semiclassically, the next assnmption that the 
constant E in Eq. (1.2), may be replaced by its time-dependent connterpart E{t), merits a 
critical examination. Strictly speaking, this replacement is invalid since as soon as the held 
acqnires time dependence, it can directly excite negative-energy particles to levels above 
the mass gap, by a pertnrbative mechanism, withont reconrse to any tnnneling or barrier 
penetration mechanism. The above replacement may be jnstihed as a kind of a local- 
held approximation provided the time-dependence of E{t) is so weak that it only excites 
freqnencies with magnitndes negligible compared to the mass gap. Eqnivalently, the held 
needs to be constant on the scale of the Compton wavelength m~^. Fields that are constant 
over snch large time intervals are nnlikely to be prodnced in heavy-ion collisions. Moreover, 
a jnstihcation for the local-held approximation, which conld be looked for a posteriori is not 
forthcoming from the explicit numerical results presented in [Q . 

The purpose of this paper is to take the time dependence of the electric held into account. 
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as properly as a semiclassical formalism would allow, and present a new source term based 
on Fig. lb, for use in the QGP transport equation, thereby avoiding the unphysical aspect 
that we have pointed out in the current formulation. The mechanism depicted in Fig. lb 


was used by Martin et ah to calculate the total energy deposited during a certain time 
interval, under certain assumptions. They did not derive a source term. 


II. NEW SOURCE TERM 

For fermion pair production in the presence of a space-time-dependent background held 
the amplitude to the lowest order (Fig. lb) is given by 

-ie u{pi) 4{K) v{p 2 ), 

and the probability for pair production is 

where 


T = tr [(A + m) 4{K) - m) 4{-K)]. 

Since our aim is to obtain the production rate in the one-particle phase space, for insertion 
in the Boltzmann equation for a one-particle distribution function f{x,p), we choose to 
perform the iF-integration hrst. Next, in order to evaluate T, we introduce 

P = Pi + P2 and p = Pi - p2, (2.2) 

and use pi ■ P 2 + xn? = P^/2. This gives 

T = (P + p) ■ A{P) (P - p) ■ A{-P) - 2P^A{P) ■ A{-P) + (P + p) ■ A{-P) (P - p) ■ A{P) 
= 2 P ■ A{P)P ■ A{-P) - 2 p2yl(P) ■ A{-P) - 2 p ■ A{P) p ■ A{-P). (2.3) 

Now recall 
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( 2 . 4 ) 


- If^{p)f‘"(-p) = P ■ MP) P ■ A-P) - pA(f) ■ A(-P) 

= |E(P)|2-|B(P)|^ 

To be able to compare our results with those based on Eq. (1.2), we ignore the magnetic 
held B(P). Note also that pair production is an electric effect |l|. Equations (2.3) and (2.4) 
give us 

T = 2 lE(P)l^-2p-A(P) p-A(-P). (2.5) 

Note that the remaining term p ■ A(P) p ■ A{—P) in Eq. (2.5) is also invariant under the 
gauge transformation 


A,{P)^A^{P)-zP^x{P), 

because p ■ P = 0. We shall now rewrite this term also in terms of E(P). To that end, we 
choose the following gauge: 

A^ = 0 = A^, 

y4"(x) = 9/3a(a;), (a,/3 = 0or3) (2.6) 

where a(x) = a(t, z) is an arbitrary real function vanishing at t and = ±cx), and is 
the totally antisymmetric Levi-Civita tensor in (1 + 1) dimensions. (This gauge choice is 
natural in the colour-plates model where E is taken along the collision or z-axis.) Equation 
(2.6) ensures the gauge condition d^A^{x) = 0. We have 

kl“(P) = (27r)"2y' d^x exp(-zP-x) ff“(x) 

= i P /3 (27r)“^ J d^x exp(—iP ■ x) a{x) 

= i P/3 a(P) 

= i P“ a(P), (2.7) 

where P is the dual of P. 
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Keeping in mind the boost-invariant central-rapidity region expected to be formed in 
ultrarelativistic heavy-ion collisions, we assume that the electric held E(a;) depends only on 
the boost-invariant variable r = \/t‘^ — z^. Hence 

E(P) = (27r)“^y d‘^x exp{—iP-x) E(a;) 

= (5(^)(Pr) J dt dz exp[-i{Pot - P,z)] E(r). (2.8) 

It follows that Pit = —P 2 T and hence the two particles have identical transverse mass, say 
rriT- The relevant components of their 4-momenta are 

= rriT{coshHi, sinh|/i), P 2 = niT{cosh 1 / 2 , sinhi/ 2 ), (2.9a) 

such that 

p" = —mT(sinhpi, coshpi), P 2 = —mT(sinhp 2 , cosh 1 / 2 ), (2.9b) 

where 1 / 1,2 are their rapidities. For a future reference we note the following results 

Pi-pi = 0 = P2-p2, Pi ■ P2 =-Pi ■ p2, (2.10) 

P-P = 0, P^ = -P^, (2.11) 

= {Pi+P 2 f = 2m\.[l + cosh(//i - Pa)]- (2.12) 

Now in order to rewrite p ■ A{P) p ■ A{—P) in terms of E(P), we use Eqs. (2.7), (2.2), 
(2.10) and (2.9) to get 

p ■ A{P) p ■ A{—P) = 4 mPp sinh^(//i — 1 ) 2 ) a{P)a{—P). (2-13) 

We also note from Eqs. (2.4), (2.7) and (2.11) that 

|E(P)p = (p2)2 a(P)a(-P). (2.14) 

Finally, eliminating a{P)a{—P) between Eqs. (2.13) and (2.14) we get the desired expression 
p ■ A{P) p ■ A{-P) = 4 sinh2(//i - 1 / 2 ) |E(P)|2 / {P'^f. (2.15) 
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The occurrence of (P^)^ in the denominator is noteworthy. It indicates a highly nonlocal 
nature of this term in the conhguration space; the production rate depends not only on the 
instantaneous strength of the held but also on its derivatives. This feature is absent in Eq. 
(1-2). Substituting P^ from Eq. (2.12) in Eq. (2.15) and then substituting the resulting 
expression in Eq. (2.5) we get 

T = 2 |E(P)|2 sech^ • 

Thus the probability W in Eq. (2.1) becomes 

In order to evaluate E(P) occurring in Eq. (2.16) we dehne (for t > 0, t^ > z^) 

T^ = t'^ — z^, 7] = tanh~^(z/t), 

M^ = P^-Pl y = tanh-i(P,/Po), 

and substitute 


dt dz = T dr dr], 

Pot — PzZ = M T cosh(y — T]) = M T cosh 6*, 

in Eq. (2.8) to get 

poo poo 

E{P) = S^'^\Pt) / d9 dr T exp^—iMr cosh 9) E{t). (2-17) 

J—oo J 0 

We shall assume for the sake of dehniteness and simplicity that E{t) depends on only one 
time scale (tq) and is of the form 

P(r) = Po exp(-r/ro), Tq > 0. (2.18) 

Self-consistent numerical calculations reported in Ref. [Q predict r-dependence which is 
consistent with Eq. (2.18). Substituting P(r) given in Eq. (2.18) in Eq. (2.17) and 
performing the r-integration by parts, we get 
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The ^-integration too can be done analytically which gives 


E(P) = E„ t| 2 / (1 + MV|), 


(2.19) 


where 


I{x) = -l-f 


In 


1 + ix + A/r+a? 


( 2 . 20 ) 


y/1 + a;2 1 + ix — VT + a? 

We continue with our derivation of W by substituting E{P) from Eq. (2.19) in Eq. 
(2.16). This entails the calculation of [5 *'^^(Pt)]^ which we perform as follows: 


[5^^)(Pt)]^ = (5^^^(Pr) J exp(iPT-xr) / (27r)^ 

^ d^XT / (27r)l 


( 2 . 21 ) 


This is a standard procedure and it allows us to write the probability per unit transverse 
area. Next, recall the occurrence of (P^)^ in Eq. (2.15) and the comment made thereafter. 
In order to recover the dependence of W on the instantaneous strength of the field |P(r)p 
and to obtain a time-dependent rate, we shall employ the integral representation: 


POO 

iro = / exp(-2r/ro) dr. 

JO 


( 2 . 22 ) 


Equations (2.19)-(2.22) together with Eq. (2.16) yield 


ToV „,2 > IP/ m2 fd^pid?p2 


dW = a ( — ) d^XT dr |P(t)| 


El E 2 


(5(2) (P^ 


X |/(Mro)|2 sech^ 


(1 + M2r2)- 


Note that 


d^pi d^p2 2 


= d Pit dyi d psT dpa- 


El E 2 

The integration over p 2 r can be done trivially with the help of the delta function = 


Pit + P2t)- In order to obtain the production rate in the one-particle phase space we follow 












the literature and assume as in Ref. that if a particle is produced at z and t, it must 
appear with the longitudinal velocity Pz/E = z/t. (This is reminiscent of the Thomas- 
Fermi approximation in condensed-matter physics.) It allows us to identify the spatial and 
momentum rapidities and we can write dr dyi = dt dz/r. 

Combining these results we hnally get the new source term: 


d^x dt d^pT 


TT 


(2.23) 

where the electric held E{t) and the function I{Mtq) are dehned in Eqs. (2.18) and (2.20), 
respectively, and from Eq. (2.12) the invariant mass of the pair, M, is given by 


M = 2 rriT cosh i(j/i — 2 / 2 )- 


III. RESULTS 

We now compare the old and the new source terms given in Eqs. (1.2) and (2.23), 
respectively. Both are invariant under Lorentz boosts in the longitudinal or 2 ; direction. In 
the limit E —>■ 0, (1.2) and all of its derivatives vanish, while (2.23) behaves like ~ In 
the large E limit, (1.2) behaves like ~ E ln(E), while (2.23) behaves like ~ Thus the 
ratio of the new and the old source terms increases with the strength of the electric held. 
The transverse-momentum {pt) dependence arises naturally in the derivation of Eq. (2.23). 
The longitudinal-momentum dependence too can be recovered, if necessary. 

We also studied these two terms numerically as a function of the proper time r and the 
transverse momentum pt, for various values of the input parameters, namely the particle 
mass m, the parameter tq appearing in Eq. (2.18), the hne-structure constant a, and the 
initial strength of the electric held Eg. Some representative results are shown in Figs. (2)- 
(4), where we have taken m = 5 MeV, tq = 0.1 fm, a = 0.2, and ^Eq = 10 GeV/fm^. 
Signihcant diherences are evident. 
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In view of the above comments and the numerical results presented here, the two source 
terms, Eqs. (1.2) and (2.23), when substituted in the Boltzmann equation are expected to 
contribute differently to the formation of the quark-gluon plasma. It is straightforward to 
extend Eq. (2.23) to the case of coloured quark and antiquark production, by introducing 
an appropriate colour factor. A similar source term for gluon production in a space-time- 
dependent colour held would also be required; that work is in progress [T^. It would be 


interesting to substitute these new source terms in a nonabelian transport equation [jl^, and 
study formation of quark-gluon plasma. We plan to undertake this work in the future. 
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FIGURES 

FIG. 1. Schematic diagrams for pair production in (a) time-independent and (b) 
time-dependent electric fields. In (a), crosses indicate interactions (to all orders) with the vac¬ 
uum and the vertical line represents the “cut”, (b) represents pair production in the lowest-order 
perturbation theory. 

FIG. 2. The old source term, Eq. (1.2), in units of fm“^ GeV“^, as a function of the transverse 
momentum pT (GeV), for various values of the proper time r. The curves are labelled by r (fm). 
Input parameters are as in the text. 

FIG. 3. Same as Fig. 2, for the new source term, Eq. (2.23). Note the difference in the scale. 

FIG. 4. The ratio of the new and the old source terms. The curves are labelled by r (fm). 
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